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RESUMEN 

El resumen sera traducido al espanol por los editores. Fast stellar winds can sweep up ambient media 
and form bubbles. The evolution of a bubble is largely controlled by the content and physical conditions of 
the shocked fast wind in its interior. This hot gas was not clearly observed until the recent advent of Chandra 
and XMM-Newton X-ray observatories. To date, diffuse X-ray emission has been unambiguously detected from 
two circumstcllar bubbles blown by WR stars, four planetary nebulae, and two supcrbubbles blown by young 
clusters. Model fits to the X-ray spectra show that the circumstellar bubbles are dominated by hot gas with low 
temperatures (< 3 x 10 6 K), while the interstellar bubbles contain significant fractions of hotter gas (> 5 x 10 6 
K). In all cases, large discrepancies in the X-ray luminosity are found between observations and conventional 
models of bubbles. Future theoretical models of bubbles need to re-examine the validity of heat conduction and 
take into account realistic microscopic processes such as mass loading from dense clumps/knots and turbulent 
mixing. Chandra ACIS-S observation of NGC 6888 will shed light on these astrophysical processes. 

ABSTRACT 

Fast stellar winds can sweep up ambient media and form bubbles. The evolution of a bubble is largely controlled 
by the content and physical conditions of the shocked fast wind in its interior. This hot gas was not clearly 
observed until the recent advent of Chandra and XMM-Newton X-ray observatories. To date, diffuse X-ray 
emission has been unambiguously detected from two circumstellar bubbles blown by WR stars, four planetary 
nebulae, and two supcrbubbles blown by young clusters. Model fits to the X-ray spectra show that the 
circumstcllar bubbles arc dominated by hot gas with low temperatures (< 3 x 10 6 K), while the interstellar 
bubbles contain significant fractions of hotter gas (> 5 x 10 6 K). In all cases, large discrepancies in the X-ray 
luminosity are found between observations and conventional models of bubbles. Future theoretical models of 
bubbles need to re-examine the validity of heat conduction and take into account realistic microscopic processes 
such as mass loading from dense clumps/knots and turbulent mixing. Chandra ACIS-S observation of NGC 6888 
will shed light on these astrophysical processes. 
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STARS: MASS LOSS — STARS: WOLF-RAYET STARS — X-RAYS 
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1. INTRODUCTION: A BRIEF HISTORY OF 
BUBBLE STUDIES 

In 1965, Johnson & Hogg reported three shell 
nebulae, NGC 2359, NGC 6888, and the nebula 
around HD 50896, and suggested that these nebu- 
lae were formed by interactions between the mass 
ejected by the central Wolf-Rayct (WR) stars and 
the ambient interstellar gas. In modern terminol- 
ogy, these shell nebulae are called "bubbles", the 
mass ejected by a WR star is "fast stellar wind" , the 
shell nebula around HD 50896 has been cataloged as 
"S 308" , and the ambient medium is designated "cir- 
cumstellar medium", i.e., ejected stellar material. 

The earliest theoretical treatments of interac- 
tions between fast stellar winds and the interstellar 
medium (ISM) were motivated by the central cav- 
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ities and high- velocity motions observed in H II re- 
gions (Mathews 1966; Pikel'ner 1968; Pikel'ner & 
Schcheglov 1969; Dyson & de Vries 1972). Mod- 
els of wind-ISM interaction were constructed specif- 
ically for shell nebulae around WR stars by Avedis- 
ova (1972), and a numerical calculation of wind-ISM 
interaction was provided by Falle (1975). 

In 1974, Jenkins & Meloy reported Copernicus 
satellite observations of 32 early-type stars, which 
revealed ubiquitous shallow interstellar O VI absorp- 
tion. To explain this O VI absorption, Castor, Mc- 
Cray, & Weaver (1975) modeled wind-ISM inter- 
action and coined the term "interstellar bubble". 
Weaver et al. (1977) elaborated on this model and 
gained popularity because of their detailed descrip- 
tion of the physical conditions and structure of a 
bubble and their specific characterization of X-ray 
emission and O VI column density expected from a 
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bubble interior. 

High-quality, high-resolution X-ray and far UV 
observations were not possible until Chandra X-ray 
Observatory, XMM-Newton X-ray Observatory, and 
Far Ultraviolet Spectroscopic Explorer (FUSE) satel- 
lites were launched. It is finally possible to observe 
the hot gas in bubble interiors and critically examine 
the bubble models. The analysis of FUSE observa- 
tions is complex, and it is pre-mature to conclude on 
the O VI absorption from bubbles. Therefore, in this 
paper we will concentrate on only X-ray observations 
of bubbles and compare them to model predictions. 

2. BUBBLE MODELS 

Most of the early models of bubbles assume that 
the stellar wind interacts with the ambient medium 
through momentum transfer and find that the shell 
expansion follows r oc t 1 / 2 , where r is the shell radius 
and t is the dynamic age (Mathews 1966; Pikel'ner 
1968; Avedisova 1972; Steigman, Strittmatter, & 
Williams 1975). Their assumption of momentum 
conservation may be convenient for 1-D calculations 
in which stellar wind docs not get deflected. How- 
ever, in the 3-D space, momentum is a vector and the 
momentum flux of a stellar wind is null (otherwise 
the star itself would experience a rocket effect and 
gain momentum); it is unphysical to assume that 
the wind momentum in every specific direction is 
conserved. 

Dyson & de Vries (1972) were the first to sug- 
gest that the expansion of a bubble is driven by the 
thermal pressure of the shocked wind in the bubble 
interior, and follows the r oc t 3 ^ 5 law. This formed 
the basis of the interstellar bubble model by Castor, 
McCray, & Weaver (1975), who implemented heat 
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Fig. 2. Temperature and density profiles of a pressure- 
driven interstellar bubble, from Weaver et al. (1977). 

conduction at the interface between the hot interior 
and the cool swept-up shell. 

The basic structure of an interstellar bubble, il- 
lustrated in Figure 1, can be divided into four zones: 
(a) freely expanding stellar wind, (b) shocked stel- 
lar wind, (c) swept-up ISM, and (d) ambient ISM. 
The temperature and density structure of a bubble 
is illustrated in Figure 2. Two shocks are present in 
a bubble, an adiabatic shock at R±, between zones 
(a) and (b), and an isothermal shock at R2, between 
zones (c) and (d). The hot, shocked stellar wind and 
the swept-up ISM are separated by a contact discon- 
tinuity at R c . In the shocked stellar wind layer, the 
temperature structure is modified by heat conduc- 
tion, and the density at the outer edge is raised by 
mass evaporation across R c from the dense swept- 
up interstellar shell. At the conduction front, the 
temperatures are high enough to produce collision- 
ally ionized +5 , which may be responsible for the 
interstellar O VI absorption detected in Copernicus 
observations of early- type stars. 

For a pressure-driven bubble in a homogeneous 
ISM, the X-ray emissivity can be integrated over the 
volume with X-ray-emitting temperatures to deter- 
mine the expected X-ray luminosity in soft energy 
band. As given by Chu et al. (1995), the expected 
X-ray luminosity in the 0.1-2.4 keV band can be 
expressed in wind parameters as Lx ~ (2 x 10 35 
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Fig. 1. Schematic structure of a pressure-driven inter- 
stellar bubble, from Weaver et al. (1977). 



£g"' " J , where £ is the metal- 
licity relative to the solar value, L37 is the mechan- 
ical luminosity of the stellar wind in units of 10 37 
ergs s _1 , no is the ambient density in cm -3 , and te 
is the dynamic age in units of 10 6 yr. The X-ray 
luminosity can also be expressed in bubble parame- 
ters as L x « (1.6 x 10 28 ergs b- 1 )^ 7 R^/ 7 V^ /r , 
where R pc is the shell radiua in units of pc, and 
V5 = 0.59i?p C /te is the shell expansion velocity in 
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TABLE 1 
X-RAY OBSERVING FACILITIES 



X-ray 


Mission 


Imaging 


Angular 


Energy 


Observatory 


Duration 


Spectrometer 51 


Resolution 


Range (keV) 


Einstein 


1978-1981 


IPC 


120" 


0.2-3.5 


ROSAT 


1990-1998 


PSPC 


30" 


0.1-2.4 


ASCA 


1993-2001 


SIS 


150" 


0.4-10 


Chandra 


1999-pres. 


ACIS 


1" 


0.1-10 


XMM-Newton 


1999-pres. 


EPIC 


10" 


0.1-15 



a IPC: Imaging Proportional Counter; PSPC: Position Sensitive Proportional Counter; SIS: Solid-State Imag- 
ing Spectrometer; ACIS: Advanced Camera for Imaging and Spectroscopy; EPIC: European Photon Imaging 
Camera 



units of km s 1 . 

Shell nebulae around WR stars are circumstel- 
lar bubbles consisting of stellar material ejected by 
the progenitors at a red supergiant (RSG) phase 
or a luminous blue variable (LBV) phase; the cir- 
cumstcllar medium is far from homogeneous. The 
formation and evolution of WR bubbles have been 
hydrodynamically simulated by Garcfa-Segura et al. 
(1996a, b) for WR stars that have evolved through 
LBV and RSG phases, respectively. These models 
follow the same basic principles as Weaver et al., 
but have encorporated realistic stellar mass loss his- 
tory in their calculations. While these models are 
successful in reproducing the nebular morphologies, 
they find that the stellar wind luminosity expected 
from the observed bubble dynamics is often more 
than an order of magnitude lower than that derived 
from observations of the wind directly. 

The clumpy morphology of WR bubbles implies 
that the stellar wind may be interacting with small 
fragments of nebular material, thus mass loading 
may be an important process that modifies the phys- 
ical conditions of a bubble interior. Adiabatic bub- 
bles with conductive evaporation and hydrodynamic 
ablation have been modeled by Pittard et al. (2001a, 
b). They find the X-ray emissivity and temperature 
profiles of a bubble interior can rise or fall toward 
the outer edge, depending on the ratio of wind mass 
to injected nebular mass. 

Wind-blown bubbles are usually associated with 
massive stars; however, the formation of planetary 
nebulae (PNe) around low- and intermediate-mass 
stars is almost identical to that of a WR bubble (e.g., 
Kwok 1983; Chu 1993). Therefore, we will include 
PNe in our discussions in this paper. 



3. X-RAY OBSERVATIONS OF BUBBLES 

X-ray observations of wind-blown bubbles and 
PNe have been made with a number of satellites, of 
which the duration of the mission, instrument used, 
angular resolution, and photon energy range are 
summarized in Table 1. Einstein IPC observations 
detected diffuse X-ray emission from only one wind- 
blown bubble, NGC6888 (Bochkarev 1988). Many 
detections of diffuse X-ray emission from interstellar 
bubbles in HII regions based on Einstein observa- 
tions, e.g., the Orion Nebula (Ku & Chanan 1979), 
the Rosette Nebula (Leahy 1985), S155 (Fabian & 
Stewart 1983), were later shown by R OSA T obser- 
vations to be spurious, as the X-ray emission was 
resolved into stellar point sources (Chu 1994). 

In Table 2, we summarize X-ray observations of 
bubbles and PNe made with "modern" X-ray satel- 
lites within the last decade. Objects with diffuse 
X-ray emission detected are listed at the top with- 
out parentheses, and the non-detections are listed at 
the bottom in parentheses. 

3.1. Diffuse X-ray Emission from WR Bubbles 

Diffuse X-ray emission has been detected from 
two WR bubbles, NGC6888 and S308. Figures 3 
and 4 show optical [O III] and ROSAT PSPC images 
of these two bubbles (Wrigge, Wendker, & Wisotzki 
1994; Wrigge 1999). Despite the limited angular res- 
olution, the PSPC X-ray image of NGC 6888 shows 
clear limb-brightening. The X-ray morphology of 
S 308 cannot be easily visualized because the shell 
rim is occulted by the circular ring of the PSPC's 
window support structure. Analyses of the spectra 
extracted from these PSPC observations show that 
the hot gas in the interiors of these two WR bub- 
bles is dominated by gas at ~1.5xl0 6 K. Wrigge 
(1999) also suggested a high-temperature, 2.8xl0 7 
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K, component in S308; however, this component 
may be contributed by the numerous unresolved 
point sources. ASCA observations of NGC6888 de- 
tected an additional plasma component at ~8xl0 6 
K (Wrigge et al. 1998). 

These observations can be greatly improved by 
Chandra and XMM-Newton. We have been awarded 
a 100 ks Chandra observation of NGC 6888 in Cy- 
cle 4 that will be made in 2003. We have obtained 
40 ks XMM-Newton observations of S308, with the 
field of view covering the northwest quadrant of the 
bubble (see Figure 5). Although these XMM obser- 
vations were affected by high background during 2 /3 
of the total exposure time, they show clearly a limb- 
brightened distribution of the X-ray emission from 
the interior of S308. Most interestingly, a spatial 
gap is present between the outer edge of diffuse X- 
ray emission and the outer edge of [O III] emission, 

TABLE 2 

X-RAY OBSERVATIONS OF BUBBLES AND 
PLANETARY NEBULAE a < b 



ROSAT 


Chandra 


XMM-Newton 




Bubbles 




NGC 6888 


Rosette Nebula 


S308 


S308 


Omega Nebula 




Omega Nebula 






(NGC 2359) 






(NGC 3199) 






(NGC 6164-5) 






(NGC 7635) 








Planetary Nebulae 




BD+30°3639 


BD+30°3639 


NGC 7009 


NGC 6543 


NGC 6543 




Abell 30 


NGC 7027 




(60+ PNe) 


(NGC 7293) 






(He 2-90) 






(M 1-16) 





a Diffuse X-ray emission is detected only from objects 
listed at the top without parentheses. 

b References for objects with diffuse X-ray emission - 
Abell 30: Chu, Chang, & Conway (1997), Chu, & Ho 
(1995); BD+30°3639: Kastner et al. (2000); NGC 6543: 
Chu et al. (2001); NGC 6888: Wrigge, Wendker, & 
Wisotzki (1994); NGC 7009: Guerrero, Gruendl, & Chu 
(2002); NGC 7027: Kastner, Vrtilek, & Soker (2001); 
Omega Nebula: Townsley et al. (2002), Dunne et al. 
(2002, in prep.); Rosette Nebula: Townsley et al. (2001); 
S308: Wrigge (1999), Chu et al. (2002, in prep.). 




Right Ascension (J2000) 



Right Ascension (J2000) 



Fig. 3 
PSPC 



Optical [O III] A5007 image (left) and ROSAT 
image (right) of NGC ( 
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Fig. 4. Optical [O III] A5007 image (left) and ROSAT 
PSPC image (right) of S 308. 
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Fig. 5. Optical [O III] A5007 image (left) and XMM- 
Newton EPIC image (right) of the northwest quadrant 
of S 308. 

and this gap might be the conduction front that has 
been long sought after! The X-ray spectrum of S 308 
(see Figure 6) is extremely soft, indicating a plasma 
temperature of only ~lxl0 6 K. 

3.2. Diffuse X-ray Emission from Planetary 
Nebulae 

More than 60 PNe have been observed by 
ROSAT, but only three nebulae show marginally 
extended X-ray emission (Guerrero, Chu, & Gru- 
endl 2000). Chandra observations unambiguously 
resolved the diffuse X-ray emission from three PNe, 
and XMM-Newton observations resolved the diffuse 
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Fig. 6. XMM-Newton EPIC-pn and Chandra ACIS-S spectra of four PNe and one WR bubble. 



X-rays from NGC7009. See Table 2 for the names 
and references of these PNe. Figure 7 shows the most 
well-resolved PN, NGC 6543 (the Cat's Eye Nebula), 
where a limb-brightened morphology is clearly seen. 

The X-ray spectra of the diffuse emission from 
four PNe are displayed in Figure 6; all show thermal 
plasma emission. The best-fits to these spectra indi- 
cate plasma temperatures of 2-3 xlO 6 K and densi- 
ties of ^100 cm -3 , from which we conclude that the 
hot gas in PN interiors is over-prcssurized and drives 
the expansion of the optical nebular shell. 
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Fig. 7. Chandra ACIS image (left) and HST WFPC2 
Ha image overlaid by X-ray contours (right) of the PN 
NGC 6543, the Cat's Eye Nebula. 



3.3. Diffuse X-ray Emission from Superbubbles 

Diffuse X-ray emission from quiescent superbub- 
bles (without recent supernova blasts) is detected for 
the first time by Chandra - the Rosette Nebula and 
the Omega Nebula (Townsley et al. 2001, 2002). The 
diffuse X-ray emission from the Omega Nebula was 



in fact detected by ROSAT, but has never been re- 
ported. Figure 8 shows the ROSAT PSPC image and 
an optical image of the Omega Nebula; the diffuse 
X-ray emission fills the entire interior of this super- 
bubble. As the ionizing cluster of the Omega Nebula 
is only ^1 x 10 6 yr old, the hot gas must be produced 
solely by fast stellar winds. 

The diffuse X-ray emission from the two inter- 
stellar superbubbles, the Rosette Nebula and the 
Omega Nebula, is qualitatively different from the dif- 
fuse emission from circumstcllar bubbles, i.e., WR 
bubbles and PNe. First of all, the diffuse emis- 
sion from the superbubbles docs not show any limb- 
brightening. Second, X-ray spectra of the super- 
bubbles show a high-temperature (~7xl0 6 K) com- 
ponent in addition to the dominant component at 
~2xl0 6 K. Finally, the densities of the hot gas in 
superbubbles are much lower than those in the inte- 
riors of WR bubbles or PNe. 



-16 00 

a 08' 

a , 
I 12 

I IB' 




18 h 22 m 00" 30" 21 "OO" 30" lO^^O" 30" 21 m 00" 30" 
Right Ascension (J2000) Right Ascension (J2O0O) 



Fig. 8. Digitized Sky Survey red-band image (left) and 
ROSAT PSPC image (right) of the Omega Nebula. 
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TABLE 3 

PHYSICAL PROPERTIES OF HOT GAS IN 
BUBBLE INTERIORS 

Bubble Type T c N c L x 





[10 6 K] 


[cm 3 ] 


[ergs 


s- 1 ] 


Planetary Nebula 


2-3 


100 


10 31 - 


-10 32 


WR Bubble 


1, 8? 


10 


10 33 - 


10 34 


Superbubblc 


2, 7 


0.1 


10 33 - 


-10 34 
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4. OBSERVATIONS VS. THEORETICAL 
EXPECTATIONS 

X-ray observations have detected diffuse emission 
from two WR bubbles, four PNe, and two young su- 
perbubbles. The typical physical properties of the 
X-ray-emitting gas derived from Chandra and XMM- 
Newton observations are listed in Table 3. The av- 
erage X-ray spectrum of a bubble is weighted to- 
ward the densest X-ray-emitting region. As shown 
in Figure 2, the density and temperature profiles in a 
bubble interior are anti-correlatcd, thus the brightest 
emission is expected to originate from regions with 
the lowest temperatures. The observed low plasma 
temperatures and limb-brightened morphology are 
fully consistent with this expectation. The observed 
variations of plasma density among the three types 
of bubbles are also qualitatively consistent with the 
expectation from mass evaporation processes, as the 
plasma density should reflect the nebular density. 

However, two outstanding discrepancies exist be- 
tween the observations and theoretical predictions 
when they are compared quantitatively. First, the 
observed X-ray luminosities are 10-100 times lower 
than those predicted by pressure-driven bubble mod- 
els using observed bubble dynamics, such as shell 
size, expansion velocity, nebular density, and stellar 
wind mechanical luminosity. Second, in a pressure- 
driven bubble model, the hot gas near the con- 
duction front is dominated by nebular mass evapo- 
rated across the contact discontinuity, but the abun- 
dances of the X-ray-emitting gas in at least two PNe 
are consistent with those of the fast stellar winds 
instead of those of the nebular shells (Chu et al. 
2001; Arnaud et al. 1996). To resolve these prob- 
lems, high-resolution observations of more bubbles 
are needed to study the detailed astrophysical pro- 
cesses, e.g., conduction evaporation, hydrodynamic 
ablation, and turbulent mixing. Our upcoming 100 
ks Chandra observation of NGC 6888 may shed light 
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